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Abstract. 

One of the primary goals when studying stellar systems with neutron stars has been to reveal the physical properties of 
progenitors and understand how neutron star spins and birth kicks are determined. Over the years a consensus understanding 
had been developed, but recently some of the basic elements of this understanding are being challenged by current observations 
of some binary systems and their theoretical interpretation. In what follows we review such recent developments and highlight 
how they are interconnected; we particularly emphasize some of the assumptions and caveats of theoretical interpretations 
and examine their validity (e.g., in connection to the unknown radial velocities of pulsars or the nuances of multi-dimensional 
statistical analysis). The emerging picture does not erase our earlier understanding; instead it broadens it as it reveals additional 
pathways for neutron star formation and evolution: neutron stars probably form at the end of both core collapse of Fe cores of 
massive stars and electron-capture supernovae of ONeMg cores of lower-mass stars; birth kicks are required to be high (well 
in excess of 100 kms~') for some neutron stars and low (< lOOkms -1 ) for others depending on the formation process; and 
spin up may occur not just through Roche-lobe overflow but also through wind accretion or phases of hypercritical accretion 
during common envelope evolution. 
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NEUTRON STAR FORMATION 
Elements of a "Standard Picture" 

Forty years after the discovery of pulsars it is fair to 
say that an immense amount of progress has been made 
in the quest to understand how neutron stars (NS) form 
and evolve in different stellar systems. Our understand- 
ing of these processes has benefited most importantly by 
the study of binary pulsars or X-ray binaries harboring 
NS with either stellar or degenerate companions (white 
dwarfs or other NS). Progress has been made both by 
studying specific observed systems as well as whole pop- 
ulations in the Galactic field or in globular clusters. These 
analyses and results have always coupled strongly to ef- 
forts in the areas of core-collapse and supernova (SN) 
explosions of massive stars, thermonuclear explosions of 
accreting white dwarfs, and nuclear evolution of massive 
stars. 

Through this diverse body of work (for which it is im- 
possible to include here all important references, but rel- 
evant overviews may be found in 0][2)) a generally ac- 
cepted understanding seemed to have emerged anchored 
on the following crucial elements most relevant to the 
discussion presented in this paper: (i) NS form at the 
end of stellar nuclear evolution as remnants of SN explo- 
sions of massive stars of ~ 1OM up to ~ 20 — 40 M ; 
(ii) NS in close binaries typically have helium-rich pro- 
genitors at the time of core collapse; these progeni- 
tors are the helium-rich cores (typically more massive 



than 2 M Q ) of original massive binary companions which 
have lost their hydrogen-rich envelopes through binary 
interactions involving (stable or unstable) mass transfer 
between binary components and/or mass loss from the 
binarjQ (iii) at the time of their formation NS acquire 
recoil kicks due to asymmetries that develop during the 
core-collapse process; typical kick magnitudes inferred 
from studies of pulsar population kinematics and the for- 
mation of binaries with NS are in the range of 100- 
300kms _1 ; (iv) NS are born with birth spin periods of 
~ 10 ms and spin down over a time span of ~ lOMyr; 
if in binaries they can spin up and reach spin periods of 
10 ms or shorter later in their life through accretion of 
material from their companion and accompanied reduc- 
tion in magnetic -field strength. 



Recent Developments and Challenges 

In recent years a growing number of NS studies have 
challenged different aspects of the "standard picture" as 
described above. Although each one of these in isolation 
may not have been powerful enough to cause a paradigm 
shift, the current set of inter-linked arguments presented 
in the literature in the past ~ 5 years has certainly at- 
tracted attention and has opened the door to potential 



In binaries with low-mass companions NS may also form through 
the accretion-induced collapse of white dwarfs, but here we focus on 
massive binaries. 



changes in our current understanding of how NS form. In 
fact, as we will discuss in what follows, these recent de- 
velopments do not actually erase or substitute our earlier 
understanding; instead they bring forward new possibili- 
ties that had not received much attention or had not even 
been raised previously. 

With a broad-brush look, two main physical processes 
lie at the core of a lot of the issues raised in connec- 
tion to NS formation and evolution: (i) electron-capture 
supernovae (ECS), and (ii) hypercritical accretion onto 
NS in common envelope phases (dynamically unstable 
mass transfer leading to common-envelope formation 
and ejection through NS inspiral towards the core of a 
massive companion). 

The possibility of ECS was originally put forward 
as a potential channel for NS formation that does not 
involve the "usual" collapse of the iron core of a massive 
star at the end of nuclear evolution |U 0. Instead 
massive ONeMg white dwarfs can reach high enough 
densities for electrons to capture onto 24 Mg and 20 Ne 
(before Ne ignition can occur (see J3]|6]|7]), causing the 
white dwarf to lose support against gravity; ultimately 
this process is expected to lead to core collapse, a SN, 
and a NS remnant. Such a collapse could occur right 
at the progenitor mass boundary between the formation 
of massive, ONeMg white dwarfs and the collapse of 
iron cores into NS: hydrogen-rich, original progenitors 
of ~ 8 — 10 M . This mass range is, however, highly 
uncertain. 

More recently Podsiadlowski et al. proposed that 
ECS events are relevant for stars of initial mass in the 
range 8—11 M in binary systems where they can avoid 
the second dredge-up phase on the asymptotic giant 
branch (AGB) (due to loss of the hydrogen-rich enve- 
lope) and form bigger helium cores; instead single stars 
of the same mass would form ONeMg white dwarfs. 
They also suggested that ECS are relatively fast explo- 
sions compared to delayed neutrino-driven explosions 
of iron cores and could naturally lead to natal kicks of 
low magnitudes. Kitaura et al. [9] modeled ECS explo- 
sions of a ~1.3M ONeMg white dwarf (formed from 
a 2.2 Mq helium core, which in turn formed from an 
8.8 M hydrogen-rich progenitor); they find that the ex- 
plosions are rather faint (low 56 Ni production), produce 
low-mass NS with baryonic mass of 1.36M Q , and are 
not prompt but still delayed and powered by neutrino 
heating. T. Janka during his talk at this conference iPTOl 
presented more recent results and commented that the 
shock is not delayed long enough for kicks in excess of 
^100 km s~ 1 to be produced. Last, Poelarends et al. ifTTI 
investigated the late evolution of solar-metallicity single 
stars on the AGB and concluded that stars in the mass 
range 9 — 9.25 Mq could experience ECS, implying an 
ECS fraction of 4—20% among all core-collapse SN, if 
one accounts for some of the model uncertainties (AGB 



stellar wind strengths and dredge-up efficiency, semi- 
convection, convective overshooting, rotational mixing; 
see |H2][T3]). In a similar study, Siess lfT3ll concluded that 
the ECS mass range is indeed very narrow for single stars 
(<1-1.5M ); at solar metallicity he finds a mass range 
of 9.5-10.5 M© which is shifted downwards by 2M if 
core overshooting is included and could also vanish alto- 
gether depending on the metallicity dependence of AGB 
winds. 

In view of all this recent work, our current understand- 
ing of ECS in connection to NS progenitors and kicks can 
be summarized as follows: they can occur in nature and 
may contribute a small but non-negligible fraction of NS 
forming with natal kicks < 100 km s~ 1 . Their progenitor 
mass range in single and binary stars is highly uncertain, 
but may be as wide as 7 — 12 M (8] [TT] Q21, and even 
reach down to 5M for low-metallicity stars H14L po- 
tentially including helium stars less massive than 2 M 
EE] [16). 

A significant part of our understanding of how NS 
form and evolve comes from analyses of binary pulsars 
with two neutron stars (NS-NS). In the "standard" NS- 
NS formation channel [ 2 1 the first-born NS experiences 
a common-envelope phase when its hydrogen-rich com- 
panion fills its Roche lobe. As a result the orbit shrinks to 
the small orbits observed, the helium-rich progenitor of 
the second NS is revealed, and the first NS has a chance 
of getting mildly recycled due to accretion in the com- 
mon envelope and/or mass transfer from the helium star 
companion. Brown fl7| challenged the survival of the 
first NS through the common-envelope phase and argued 
that the NS will collapse instead into a black hole due to 
hyper-critical accretion (exceeding the photon Edding- 
ton limit by many orders of magnitude). He suggested 
an alternative formation scenario that envisions evolu- 
tion through a double-common-envelope phase reveal- 
ing two helium-rich NS progenitors in a tight orbit, fol- 
lowed by two consequent SN explosions. This channel 
requires that the two hydrogen-rich massive progenitors 
have almost equal masses (within 4%) and could also po- 
tentially explain the almost equal NS masses in NS-NS. 
We note however, that hydrodynamical studies of hyper- 
critical accretion (e.g., Ifl8l ) and population studies of 
NS-NS formation that have accounted for hypercritical 
accretion (e.g., [20]) have concluded that if the maximum 
NS mass is ~ 2 M or higher, then a significant fraction 
of NS can still survive the common-envelope phase and 
avoid collapse into a black hole. 

PROGENITOR MASSES 

For almost two decades now, it has been widely accepted 
that there is a minimum mass a helium-rich star must 
have for the formation of an iron core and NS formation 



at core collapse. This minimum value is somewhat de- 
pendent on the detailed treatment of convection and late 
phases of nuclear burning in stellar models, but has been 
known to be in excess of 2M (e.g., in the range 2.1- 
2.3 M , see Habets ET\ ). 

One pathway for helium stars to be less massive than 
the above limit at the time of core collapse is provided by 
mass transfer onto its binary companion, which however 
does not abort iron-core formation. Recently, several in- 
dependent groups have analysed the late stages of evolu- 
tion of NS-NS progenitors and all reached the same con- 
clusion: at least two of the known close NS-NS binaries 
(PSR J0737-3039 and PSR B 1534+12) were experienc- 
ing Roche-lobe overflow mass transfer from the helium- 
rich progenitor of the second NS onto the first NS at the 
time of the second SN explosion Il22l l23l l24l l25l l26l . 
Although the mass accretion onto the NS may be rather 
small (and yet still adequate for mild recycling), the 
mass of the helium star can be reduced significantly 
(~ 0.5 — 1.5 M from its original mass by the time of 
its core collapse E71 |28l |29l . 

Alternatively (or even concurrently?) NS formed in 
ECS are also expected to originate from low-mass he- 
lium cores. This is primarily because of the connection 
to an ONeMg core that just barely exceeds the Chan- 
drasekhar mass when electron capture and loss of sup- 
port ensue. Although the ECS progenitor mass range is 
highly uncertain (see previous section), it tends to lie be- 
low 1 1-12 M and potentially includes helium-rich pro- 
genitors less massive than 2 M (see ECS discussion in 
|3Q||3T 'I. Such low-mass progenitors would also naturally 
produce low-mass NS with baryonic masses very close 
to the Chandrasekhar limit and gravitational masses of 
~ 1.25M® (UED. 

The discovery of the double pulsar l32l and the 
mass measurement of pulsar B (the second-born NS) at 
1.25 M led Podsiadlowski et al. 11331 to suggest an ECS 
origin. A number of independent groups analyzed the 
evolutionary history of this system, given the observed 
characteristics, and derived constraints on the immedi- 
ate (right before the SN) progenitor of pulsar B. Piran 
& Shaviv Il34l were the first to allow as part of the pri- 
ors in the statistical analysis helium star progenitors less 
massive than 2 M ; they found that progenitor masses in 
the range 1.25-1.75 M were favored by their analysis 
at the 30% confidence level, while masses up to 2M 
could just be reached at 97% confidence level (see their 
Figure 2). Willems et al. [35] presented a detailed anal- 
ysis that explored the influence of all uncertainties asso- 
ciated with the systems (most notably its age and radial 
velocity with respect to the Sun); when allowing for low- 
mass progenitors they too concluded that the most likely 
helium star masses lie in the range 1.35-1.95 M , while 
the 90% confidence intervals extended up to 2.5 M (see 
their Figure 6 and Table IV). More recently, Stairs et al. 



11261 addressed some of these questions using updated 
proper motion measurements, and they too concluded 
that the immediate progenitor of pulsar B was likely less 
massive than 2 M with most likely values at ~ 1.5 M 
(see their Figure 1). These conclusions are very different 
from those of a very similar analysis of the NS-NS sys- 
tem PSR B 1534+12 for which Stairs et al. (26) report a 
most likely value of 2.45 M and a 95% confidence in- 
terval 1.6-3.9 Mq. 

Based on these results, it appears most probable that 
the immediate progenitor of pulsar B was less massive 
than 2M (most likely ~ 1.5 M ), even though this is 
clearly not universal among close NS-NS binaries. The 
next question becomes: what does this low mass tell us 
about NS formation? We would argue that the answer is 
still unclear. It is indeed very tempting to associate this 
mass with an ECS event and declare the case closed, and 
this could very well be the answer. However, we should 
not eliminate from our thinking that PSR J0737-3039 has 
also been found by several independent studies to be in 
a phase of Roche-lobe at the time of pulsar B formation. 
Consequently, we cannot exclude that the low progeni- 
tor mass is due to mass transfer (albeit with small accre- 
tion onto pulsar A) initiated before the SN explosion. We 
conclude that both of these possibilities remain as poten- 
tially relevant to the double pulsar formation. 

At present no obvious secure observational test to dis- 
tinguish between these two possibilities appears to exist. 
However we suggest that two different types of theoret- 
ical studies could shed some light and allow us to make 
further progress: (i) detailed mass transfer sequences that 
can explain not just the low progenitor mass but also the 
orbital characteristics right at the explosion time as well 
as the measured NS masses; (ii) population studies sim- 
ulating evolution forward in time focusing on the ques- 
tion of the relative frequency of ECS vs. mass transfer 
effects depending on assumptions about the highly un- 
certain ECS progenitor mass range. 

SUPERNOVA KICK MAGNITUDES 

For over more than three decades now, estimates of kick 
velocity magnitudes imparted to NS at birth have been 
obtained through the analysis of the kinematics of sin- 
gle pulsar samples of varying sizes and reliability of 
proper motion measurements (starting with Gunn & Os- 
triker [36]). Arzoumanian et al. [37| analyzed a sample 
of 79 single pulsar proper motions and concluded that 
the NS kick magnitude distribution is best described by 
a two-component Maxwellian velocity distribution with 
characteristic velocities (three-dimensional velocity dis- 
persions) of90kms~ 1 and500kms _1 with almost equal 
relative weight. More recently, Hobbs et al. [38] adopted 
an updated electron density model and used a sample of 



73 proper motion measurements for pulsars with short 
characteristic ages (< 3Myr); they concluded that the 
sample is best described by a single Maxwellian with 
a one-dimensional velocity dispersion of 265 km s , 
shifting the mean kick velocity to values ~ 400 km s . 
Very similar results were also obtained by Zou et al. |f39l. 
It is recognized that these results could be somewhat con- 
taminated by pulsars originating from disrupted binaries 
and the contribution of orbital velocities, and by a poten- 
tial bias against low velocities due to the associated small 
and hence not easily measurable proper motions. Never- 
theless, it is clear that overall these analyses point to sig- 
nificant kick magnitudes of a few hundreds of km s~ 1 . 

Additional constraints on NS kick magnitudes have 
been obtained by the analysis of the evolutionary his- 
tory of specific binary systems with NS known in the 
Galaxy, most notably binary pulsars and X-ray binaries. 
Although the majority of such studies in the past reached 
conclusions in support of significant NS kicks, more re- 
cently a number of arguments and analyses point towards 
a need for at least some NS receiving low-magnitude 
kicks. In what follows we summarize such evidence and 
place it in the context of our developing current under- 
standing. 

Neutron Stars in Globular Clusters. Perhaps one of 
the earliest arguments in favor of modest or negligible 
NS birth kicks has been anchored on the need for NS 
retention in globular clusters. The observed low-mass X- 
ray binaries (persistent and transient) and binary pulsars 
in Galactic globular clusters have led to estimates for the 
total number of NS currently present in clusters in the 
range 300— 1,000 ll40ll . implying retention fractions of 
1 — 10%. Given that current estimates of cluster escape 
velocities lie in the range ~ 20 — 50 km s~ 1 , it is evident 
that clusters could indeed pose a requirement for low 
kick magnitudes. Pfahl et al. |41 1 concluded that reten- 
tion fractions of 1-10% can be achieved if kick magni- 
tude distributions are Maxwellians with one-dimensional 
dispersions smaller than ~ 200kms~ 1 . They also found 
that the vast majority of NS are retained in massive bi- 
naries in the clusters. However the increased character- 
istic kick magnitudes derived by Hobbs et al. Il38ll re- 
emphasized the issue of NS retention in clusters. Most 
recently, Ivanova et al. lf3Tl l42l and Kuranov & Post- 
nov ll43ll concluded that the required retention fractions 
can still be obtained, if a significant fraction of NS form 
through ECS events, and if these NS acquire either no 
kicks or kicks reduced by a factor of 10 relative to kicks 
imparted to NS formed through iron core collapse. 

High Mass X-ray Binaries. Pfahl et al. 144] high- 
lighted the existence of a subclass of high-mass X-ray 
binaries characterized by relatively low orbital eccentric- 
ities (< 0.2) and wide orbits (periods in excess of 30 d). 
They argued that their formation appears highly unlikely, 
unless some NS acquire only small kicks of ~ 20 km s -1 



(see their Figures 3 and 5). Subsequently Podsiadlowski 
et al. ID associated these low kicks with NS formation 
in ECS events and accounted for the subclass of low- 
eccentricity, wide high-mass X-ray binaries assuming a 
rather wide ECS mass range of 8 — 14 M . 

Double Neutron Stars. The pattern of observed sys- 
tems requiring both low and high NS kick magnitudes 
appears to be even more perplexed for the case of NS- 
NS binaries. To our knowledge van den Heuvel ll45l was 
the first to point to NS-NS and argue in favor of low 
NS kicks based on the occurrence of small eccentrici- 
ties (< 0.3) for 5 of the 7 systems known in the Galactic 
field. Following Pfahl et al. [41] and Podsiadlowski et al. 
(HI, he suggested that some second-born NS must have 
received very small kicks (~ 20kms~ 1 ) because of ECS 
formation. On the other hand, Chaurasia & Bailes RSI 
and Ihm et al. Il47l showed that orbital evolution through 
gravitational radiation alters the eccentricity distribution: 
many of the high-eccentricity systems are drained from 
the population because they evolve fast and either merge 
and/or shift to lower eccentricities through circulariza- 
tion. 

McLaughlin et al. [48 1 highlighted a correlation be- 
tween the NS-NS eccentricities and the spin periods of 
the first-born NS in the systems, while Faulkner et al. 
Il49l argued that such a correlation would naturally arise 
if all second-born NS were formed symmetrically and no 
kicks were imparted to them. Following this initial sug- 
gestion, Dewi et al. [50| performed population studies 
adopting a simple spin-up model for the first-born NS 
and concluded that such a correlation appears only if all 
second-born NS acquire kicks from a Maxwellian with 
a velocity dispersion of < 50kms . Willems et al. lIBTI 
explored the potential role of polar kicks in forming this 
correlation and similarly conclude that low kicks for the 
second-born NS seem to be needed in order to explain 
this correlation. 

Not surprisingly, the particular case of the double pul- 
sar (PSR J0737-3039) has attracted a lot of attention in 
this discussion of NS kick magnitudes. The binary pul- 
sar has a low eccentricity (0.09) and a low second-born 
NS mass (1.25 M ), and most recently a low proper mo- 
tion (lOkrns -1 ) was reported [52]. Hence, it appears as 
the prime candidate for a system where the second-born 
NS was formed through an ECS event and acquired a 
small kick. Willems et al. 11241 assumed the "standard" 
minimum progenitor mass of 2M and adopted a flat 
prior distribution for the unknown radial velocity and 
found most probably kick magnitudes in the range 100- 
300kms _1 . Podsiadlowski et al. (331 used the low mass 
of pulsar B as the "tell-tale sign" for ECS and hence a 
small kick imparted to it, anticipating a small proper mo- 
tion measurement. Piran & Shaviv [34] used the mea- 
sured orbital characteristics and an earlier estimate of 
the proper motion to derive a probability distribution for 



the birth kick to pulsar B. In their calculation they al- 
lowed for the possibility of very low-mass progenitors 
and also made a number of restrictive assumptions about 
the system's age and its Galactic motion (the latter al- 
lowed them to by-pass the question of the unknown ra- 
dial velocity). They primarily highlighted their result de- 
rived at 30% confidence level favoring kick magnitudes 
either < 30kms _1 or in the range 50 — 80kms~'. At a 
safer confidence level of 97%, their results set an upper 
limit of 125 km s _1 (see their Figure 2C). 

Following these studies, Willems et al. [35 1 and Stairs 
et al. [26 1 presented the two most comprehensive analy- 
ses of this system to date (the latter of the two using the 
most recent and reliable measurement from pulsar timing 
of the proper motion by Kramer et al. [52]). Both reached 
rather similar conclusions for the most likely kick magni- 
tude imparted to pulsar B: 50-170 km s _1 (Willems et al. 
11351 ) and 40-80 km s" 1 (Stairs et al. EH). As expected 
these magnitudes reach much higher values if wider 
confidence intervals are considered. Willems et al. 051 
found up to 300kms _1 at 90% confidence level, assum- 
ing a Gaussian distribution for the unknown radial ve- 
locity with dispersion of 130kms _1 (motivated by pop- 
ulation synthesis calculations). Stairs et al. [26] found 
kicks in the range of 45-1000 km s" 1 for the same radial 
velocity assumptions, and of 20-140 km s , assuming 
velocity vectors are isotropically distributed around the 
transverse velocity. The issue of the radial velocities has 
been identified by these two studies in particular as quite 
important for the best possible determination of the kick 
magnitude. To avoid a long digression we discuss this is- 
sue in more detail in the Appendix where we show the 
simple assumption of isotropicity for the system velocity 
to not be justified. 

It is interesting to compare the Willems et al. 05l 
and Stairs et al. 12611 results to those obtained by Pi- 
ran & Shaviv 1341 . Although each study has used some- 
what different observational constraints as they became 
available, one can say that at zeroth order all three stud- 
ies agree that, once immediate progenitor masses lower 
than 2M are allowed as part of the priors in the sta- 
tistical analysis (motivated by ECS events or collapse 
of rather small iron cores), then kick magnitudes in the 
range 10— lOOkms -1 are clearly favored at low confi- 
dence intervals (~ 30%). For these magnitudes, Willems 
et al. Il24l and Stairs et al. [26 1 predicted a spin misalign- 
ment angle for pulsar A smaller than 20° in very good 
agreement with the upper limit of 10-15° at similar con- 
fidence intervals just recently reported by Ferdman et al. 

Some subtleties however do arise at a more detailed 
level of examination. Although these do not alter the 
basic conclusions in any way, it is somewhat important 
to briefly identify some of them, to avoid propagation 
to misunderstandings or misconceptions. For example: 



(i) a NS progenitor of mass lower than 2M does not 
necessarily imply that kicks must also be very small 
(~ 20kms~ 1 ); (ii) the current position of the system 
close to the Galactic plane does not necessarily imply 
that the system and kick velocity must be small (see the 
point made by Figure 9b in (351 ); (iii) the derived prob- 
ability distributions are inherently multi-dimensional (5- 
D: progenitor mass, orbital separation, kick magnitude, 
and two kick direction angles) and statements about pre- 
ferred kick magnitudes are quite sensitive to projection 
effects in different dimensions because of the complex 
structure of the multi-dimensional surface (as shown in 
Table III in (35 j, the statement of "most likely" kick can 
vary by a factor of ~ 7 from 10 to 75kms~' depending 
on the projection dimensionality); (iv) the system age is 
one of the poorest constrained priors in the determination 
of the probability distributions (as seen from Table III in 
Willems et al. (35), the "most likely" kick can vary by a 
factor of ~ 30 from 5 to 170 km s~ 1 depending on the as- 
sumed system age); (v) the most likely birth location of 
the system can be several kpc away from its current po- 
sition (see Figure 8 in ll35l ). so that following the three 
dimensional motion of the binary is important for the dy- 
namics. 

The results obtained for the other NS-NS binaries are 
particularly interesting in comparison to the case of the 
double pulsar. Just like with the second-born NS pro- 
genitor masses, PSR B1534+12 paints a rather differ- 
ent picture for NS kicks. Independent studies Il24l l25l 
have concluded that a minimum kick of ~ 200 km s -1 
is needed to explain its properties along with its signif- 
icant spin-orbit misalignment. It is further interesting to 
point out that this system has a rather small eccentricity 
(<0.3) and exhibits rather well the correlation between 
pulsar spin period and orbital eccentricity. Both of these 
attributes have been associated with second-born NS get- 
ting no or minimal (~ 20kms~') kicks B31 l50l fSTl . 
PSR 1534+12 therefore clearly presents a contradictory 
example. PSR B1913+16 also clearly requires a signifi- 
cant kick magnitude of at least 260 km s _1 ll24l l54l and 
references therein]. 



SPIN UP OPPORTUNITIES FOR 
DOUBLE NEUTRON STARS 

It has been long recognized that first-born NS in NS-NS 
must have been spun-up (recycled) at some stage prior to 
the second NS formation event. However, the measured 
spin periods are all rather long in the range 10-100 ms, 
in excess of the "classical" millisecond pulsars, and these 
NS are thought to be only mildly recycled. Nevertheless, 
even mild recycling requires some level of accretion onto 
the NS at some stage between the two SN events. This 



requirement is sometimes ignored in discussions of NS- 
NS formation, but can provide us with some insights and 
constraints. 

In the "standard" NS-NS formation channel, the first- 
born NS has in principle two opportunities for recycling: 
either during the common-envelope inspiral involving 
the hydrogen-rich progenitor of the second NS or dur- 
ing a mass-transfer phase driven by Roche-lobe over- 
flow from the helium core exposed after the common- 
envelope phase. As mentioned already, the progenitors 
of at least two of the known NS-NS binaries were indeed 
in such a mass transfer phase from the helium-rich com- 
panion to the first-born NS at the time of the second SN 
explosion. Therefore the need for NS spin-up during the 
common envelope is alleviated for these systems. 

It is also interesting that all explanations of the corre- 
lation between spin period and eccentricity rely on a (lin- 
ear) scaling between the accretion amount responsible 
for the spin-up and the mass of the second-born NS pro- 
genitor right before the explosion. Such a scaling arises 
naturally if spin-up happens during a phase of Roche- 
lobe overflow mass transfer from a helium-rich compan- 
ion and is supported by mass-transfer calculations of NS 
helium-star binaries [27 28, 29] (see also Fig. 2 in [51]). 
Willems et al. [15 1] furthermore noted that the presence of 
high-eccentricity systems on the spin period-eccentricity 
relation could support the existence of a similar scaling 
for NS accretion during common envelope phases. 

Last, it is important to realize that in the case of the 
double common-envelope channel put forward by Brown 
fPTl . the opportunities for NS spin-up are limited and po- 
tentially negligible. Given how similar the masses of the 
two progenitors are, the time between the two SN ex- 
plosions is very short (< l,000yr) and it is difficult to 
imagine that in this short time both the orbital and stellar 
characteristics would happen to be appropriate for estab- 
lishing Roche-lobe overflow and spinning up the first- 
born NS. In their population simulations Belczynski et 
al. EOl did indeed identify this NS-NS formation chan- 
nel, but, for a maximum NS mass of 3 M , it only con- 
tributed 8% of the total NS-NS population. Most impor- 
tantly though Belczynski et al. did not find any oppor- 
tunity for NS recycling. They also considered a model 
with a maximum NS mass of 1.5 M fi (model D2) and 
naturally found a significant contribution to NS-NS for- 
mation from the double common-envelope channel, but 
the issue of recycling remained. More recently Dewi et 
al. Il55l modeled NS-NS formation through this channel 
and obtained rates very similar to those by Belczynski et 
al. (20], but did not comment on the difficulty with recy- 
cling the first-born NS in the systems. 

At present we would say that the question of when and 
under what circumstances first-born NS in NS-NS bina- 
ries get recycled is rather unclear. If the double common- 
envelope channel is shown to be dominant in some way, 



then spin-up during mass transfer from the helium-rich 
companion would be necessary. It is possible that accre- 
tion during common envelope inspiral or wind accretion 
alone without Roche-lobe overflow would be sufficient to 
recycle the first born NS, although so far neither of these 
have been examined quantitatively. Also, if the correla- 
tion between spin period and orbital eccentricity persists, 
the amount of mass accreted by the NS must scale with 
the donor mass regardless of the mode of mass transfer. 

CURRENT CONCLUSIONS: ONE SIZE 
DOES NOT FIT ALL 

The picture of NS formation from 10 years ago has 
clearly changed and paradigms have shifted - NS forma- 
tion now appears more diverse than we thought both in 
terms of progenitor masses and birth kicks, and thus also 
in terms of core-collapse physics, asymmetry develop- 
ments, and maybe even kick mechanisms. 

In the last 3 years, the number of papers favoring small 
or even zero NS kicks has mushroomed. A lot of this re- 
cent work is driven by the increasing observational in- 
formation regarding binary pulsars as well as new de- 
tections. We summarize the current picture as follows: 

(1) NS form from the iron core collapse of helium-rich 
progenitors with masses above ~ 2M , though lower 
masses may be possible if significant mass loss through 
Roche lobe overflow occurs shortly before NS formation. 
Alternatively, NS form from lower-mass helium stars 
through ECS of massive ONeMg cores. The ECS pro- 
genitor mass range is currently highly uncertain, but the 
channel may be particularly accessible for NS formation 
in binaries where the NS progenitors lose their hydrogen- 
rich envelope before the second dredge-up on the AGB. 

(2) Typical NS kicks of 200-300 km s _1 are still favored 
by analyses of pulsar kinematics and are required for 
many of the known binary NS systems. However sub- 
populations of or some specific NS binaries clearly fa- 
vor smaller kick magnitudes (lower than ~ 100 km s -1 , 
although probably not as low as ~ lOkrns -1 (T. Janka 
2007, private communication); these may be naturally 
provided during ECS events because of the short duration 
of the explosion limiting the development of asymme- 
tries and build-up of larger kicks. (3) NS recycling back 
to ~ms spin periods occurs through long-term, steady 
accretion through Roche-lobe overflow. For some of the 
known double NS systems such mass transfer from the 
helium-rich progenitor of the second NS is actually re- 
quired; however, depending on the dominant formation 
channel(s), mild recycling to ~ 10 ms may be possible 
(or even required) during the common envelope phase or 
through wind accretion from the helium-rich NS com- 
panion. The diversity is apparent and it is providing us 



with new clues for the formation and origin of NS (see 
also 041 for a related discussion). 

In the near future we look forward to more analyses 
that can hopefully constrain some of the current uncer- 
tainties: what is the progenitor mass range for NS for- 
mation through both iron-core collapse and ECS, espe- 
cially in binaries over a range of metallicities? Is there a 
minimum and/or maximum NS kick associated with ECS 
events? Do NS survive hypercritical accretion in com- 
mon envelopes? Which progenitor properties determine 
the birth NS masses? 
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APPENDIX 

All of the analyses of evolutionary histories and kine- 
matics of pulsar systems suffer from the lack of any con- 
straints on their radial velocities 11241 l25l l56l [35l |26|| . In 
the case of PSR J0737-3039, probably the most exhaus- 
tively studied binary pulsar, the question of the appro- 
priate radial velocity distribution that should be adopted 
is one of the most important remaining uncertainties in 
the models (along with the age). Piran & Shaviv 0411571 
did not assume anything, but therefore had to restrict 
the Galactic motion to one dimension; Willems et al. 
E4l assumed a broad flat radial velocity distribution; 
Willems et al. 051 used population synthesis for the for- 
mation and motion of NS-NS binaries to derive theoreti- 
cal radial velocity distributions filling in the void left by 
the lack of empirical constraints; Stairs et al. Il26l also 
adopted the latter in combination with updated proper 
motion constraints, but at the same time argued in fa- 
vor of a radial velocity distribution obtained by assum- 
ing an isotropic distribution for the angle 9 between the 
full 3-D space velocity and our line-of-sight. The choice 
between these two priors alters their most likely kick ve- 
locity at 95% confidence level from 45-1005 km s -1 to 
20-140 km s . Given this level of influence on results, 
we examine this question in more detail in what follows. 

The isotropic assumption appears reasonable in the ab- 
sence of other constraints and most importantly appears 
intuitive. Quantitatively, it leads to a much narrower ra- 
dial velocity PDF than the one obtained from population 
synthesis models (see Figure 1 in 11261 ). However, it is 



possible that the motion of the NS-NS population intro- 
duces velocity anisotropicities, especially when proper 
motion and location constraints are incorporated. In par- 
ticular, if the motion is dominated by Galactic rotation, 
as it is for the progenitors of NS-NS, the systemic ve- 
locity distribution will be strongly anisotropic and de- 
pend on the direction of the line-of-sight. Specifically 
let us consider the following thought experiment: if only 
small kicks are imparted to the second-born NS at birth, 
an observer located at the Galactic center will see a 9- 
distribution that is strongly peaked near angles of 90° . It 
is hard to envision that such a distribution would appear 
isotropic when the observer's viewpoint is shifted from 
the Galactic center to the position of the Sun. 

To illustrate this further, we use the theoretical NS- 
NS populations synthesized by Willems et al. 051 in 
the case of isotropic SN kick velocities with magnitudes 
drawn from Maxwellian distributions with dispersions of 
50km s _1 (as an example of pre-dominantly low kicks) 
and 260kms _1 (as inferred by Hobbs et al. P8l ). We 
follow the formation and motion of these systems in 
the Galaxy and use the results to derive present-day 
radial and transverse velocities of NS-NS binaries as 
well as PDFs for the angle 9 between the 3-D space 
velocities and the line-of-sight. We note that with this 
definition for 9 the radial velocity V,- is related to the 
transverse velocity V t as V r — V f /tan0, so 9 = points 
away from us and 9 = n towards us. Stairs et al. Il26l 
used this relation to derive an empirical radial velocity 
distribution for PSR J0737-3039 using V, = lOkms" 1 
1 52 1 and assuming an isotropic distribution for 9. 

In Figure [T[ we show the resulting present-day PDFs 
of the angle 9 for the entire NS-NS population in the 
two kick velocity models, for the sub-population of sys- 
tems within 50 pc from the Galactic plane, and for the 
sub-population of systems within 50 pc from the Galac- 
tic plane with a proper motion of less thanlOkms -1 . 
The distributions are clearly anisotropic for both kick 
velocity models, with a strong preference for directions 
along the line-of-sight (towards or away from us) emerg- 
ing when the population of NS-NS systems is restricted 
to binaries within 50 pc from the Galactic plane and a 
small proper motion of less than lOkrns -1 , characteris- 
tic of PSR J0737-3039. In Figure|2] we show the 0-PDFs 
for the sub-population of NS-NS binaries located within 
±20° from PSR J0737-3039's position in Galactic longi- 
tude and latitude. Here too, the distribution of 9 -angles 
is clearly anisotropic, with velocity vectors pointing pri- 
marily towards us. The anisotropy persists if the sam- 
ple population is restricted to systems within 1 kpc from 
the Sun (the current upper limit on the distance from the 
Sun to PSR J0737-3039 ||52]). The same holds true when 
combining the constraints on the vertical distance to the 
Galactic plane, the proper motion, and the position in 
terms of Galactic coordinates, though the statistics of the 




FIGURE 1. Present-day distributions of the direction angle 6 between NS-NS systemic velocities and our line-of-sight, for 
population synthesis models adopting Maxwellian kick velocity distributions with dispersions of 50kms _1 (left) and 260kms~' 
(right). Solid lines represent the full population of NS-NS binaries, dashed lines the sub-population of systems within 50 pc from 
the Galactic plane (CutO), and dot-dashed lines the sub-population of systems within 50 pc from the Galactic plane with a proper 
motion of less than 10 km s (Cutl). For comparison, the isotropic 6 -distribution f{6) = (sin0)/2 is also shown (dotted line). All 
PDFs are constructed using the NS-NS populations synthesized by Willems et al. (2006). 
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FIGURE 2. As Figure [7] but with different constraints imposed on the population. Solid lines represent the full population of 
NS-NS binaries (as in Figure [TJ, and dashed lines the sub-population of systems located within ±20° from PSR J0737-3039's 
position in Galactic longitude and latitude (/ = 245.2°, b = —4.5°; Cut3). The dotted line represents the isotropic 0-distribution 

/(0) = (sine)/2. 



simulated samples become considerably poorer. 

We conclude that, although it appears intuitive that 
the direction of the space velocity of the double pulsar 
is distributed isotropically around the line-of-sight, the 
system's motion due to Galactic rotation and the post- 
SN peculiar velocity combined with the non-central posi- 
tion of observers on Earth induce significant anisotropics 
in the double pulsar's most likely space velocity. Such 
anisotropies are implicitly taken into account in the theo- 
retical radial velocity distributions derived and used by 
Willems et al. [35] to constrain the formation of the 
second-born pulsar in PSR J0737-3039; therefore the re- 
sulting broader distributions of radial velocities are phys- 
ically motivated and self-consistent compared to the nar- 
rower distribution obtained with the ad hoc assumption 
of isotropic velocities with respect to our line-of-sight. 
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